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(g) Time domain filtering for NMR phased array Imaging. 



I A method and apparatus for combining NMR response data of a sample from a plurality of closely 
spaced RF receiver colls of an NMR phased array in the tjme domain to fomn a composite NMR image 
wherein each of the RF receiver coils receives a different respective one of a plurality of NMR response 
signals, each of which is evoked from a portion of the sample within a field of view of a respective one of 
the receiver coils. The response signals are conditioned to develop a plurality of data point signals 
corresponding to the magnitude of each of the respective response signals from each of the receiver 
coils at successive time intervals. The data point signals are convolved by a time domain representation 
of a field map of the respective one of the receh/er coils generating the corresponding one of the 
response signals. The convolved signals are combined on a time domain point-by-point t>asis to 
produce a time domain representation of tiie composite NMR image of tiie sample. 
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BACKGROUND OF THE INVENTION 

The present invention relates to nuclear magn tic resonance (NMR) imaging and, more particularly, to 
methods and apparatus for combining the simultaneously received data from a plurality of radio-frequency (RF) 

5 coils of an NMR phased array in the time, rather than image, domain to produce a composite image having 
high signal-to-noise ratio (SNR) throughout the image. 

The tenn "NMR phased anray" refers to apparatus, such as shown in Roemer et al. U.S. Patent No. 
4,871,969 (the disclosure of which is incorporated herein by reference), wherein a plurality of closely-spaced 
RF coils is employed for simultaneously receiving different NMR response signals from associated portions of 

10 a sample (such as a patient in medical imaging) and combining the separate data from each coil to produce a 
single composite NMR image of the sample. By overlapping adjacent coils and connecting each coil to the input 
of an associated low^input-impedance preamplifier channel, the high SNR of a single surface coil can be main- 
tained over fields-of-view (FOV) characteristic of remote coils. 

Currently, composite images for NMR phased arrays are reconstructed in the innage domain by combining 

15 the individual image contributions on a weighted, point-by-point basis after first acquiring the complete NMR 
images for each separate coil. The reason for acquiring the separate images firet is that the optimum set of 
weights needed to maximize SNR when combining the separate signals to produce the composite image is a 
function of position, and so varies ttom point to point While the phase shifters and transformers of the setup 
shown in Rgure 6 of Roemer et al. U.S. Patent No. 4.826,162 can be adjusted to provide a composite imag 

20 in the time domain having a high SNR at any particular point, different weighting must be applied for each point 
in order to obtain good sensitivity over the whole image . Thus, the conventional approach ts to first separately 
acquire the different NMR image from each coil before combining the different Individual Images, on a point- 
by-point basis, to form the composite image. 

NMR phased array imaging as described in the '1 62 patent, therefore, has the drawbacks of requiring large 

25 amounts of memory to store the separate coa Images before reconstmction and of necessitating long time 
delays between acquisition of the last data point and onset of the firet display of the reconstructed image. 

It is desirable in NMR phased array imaging to be able to combine the data from the separate receiver coils 
as it is acquired on a time domain, rather than image domain, basis without sacrificing SNR resolution. Com- 
bining the data as acquired will reduce the total memory requirements of the system since only one combined 

30 data set would have to be stored and, because only the combined data set will have to be transformed at the 
end of scanning, will also reduce the time between end-of-scan and firet appearance of the composite image. 

BRIEF SUMMARY OF THE INVENTION 

35 Among the several objects of the present invention will be noted the provision of a method and apparatus 

for forming a composite NMR image with high SNR throughout the image; the provision of a method and 
apparatus for NMR spectroscopy and NMR imaging using data combination in real time; and the provision of 
a method and apparatus which overcomes the aforementioned and other disadvantages of the prior art. 

In accordance with the invention, a method is provided for combining the simultaneously received different 

40 NMR response signals from a plurality of closely-spaced, overiapping RF receiver coils of an NMR phased array 
in the time domain, to form a composite image that has high SNR throughout the image. A filter scheme is 
utilized to develop a composite data set in the time domain, wherein each time point of the composite data is 
fonned on the basis of contributions from previous data points and future data points. The data is passed 
through filter arrangements having one-, two- and three-dimensional flltere before the signals are summed 

45 together. Each filter dimension corresponds to filtering in one of the time dimensions of k-space. i.e.. the readout 
direction; the phase encode direction; and. in the case of three-dimensional imaging, the second phase encode 
direction. RIter coefficients are chosen to combine the data in a way that is simultaneously optimal for providing 
a high SNR at multiple points of the composite image. With more tenrts added to the filter, the SNR can be 
optimized over the entire Image. 

50 In a preferred embodiment, described in greater detail below, the filter functions are determined from th 

RF magnetic field profiles of the receiving coils. 

BRIEF DESCRIPTION OF THE DRAWINGS 

55 Figure 1 (prior art) is a schematic view of an anrangement employed in th conventional method to combin 

the signals from the overiapping coils of an NMR phased array using image domain data processing techniques. 
Figure 2 is a field map sensitivity profile of a representative surface coil in the arrangement of Figure 1. 
Figure 3 is a magnitude contour plot of the time domain or k-space representation of the field map of Figure 
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Figure 4 is a filter function derlv d f r the coil of Figure 2 after Fourier transformation into Image space of 
the truncated representation of Figure 3. 

Figure 5 is a schematic view of an arrangement employed in a method of combining the coll signals of an 
NMR phased array in the time domain using filters in accordance with the present Invention. 

Figure 6 is a pictorial representation of truncated convolution in the arrangement of Figure 5 for a two- 
dimensional single slice image. 

Figure 7 is a view showing a direct implementation of the convolution sum given by equation (6), below. 

Figure 8 shows an alternative emt>odiment of the implementation shown In Figure 7. 

Figure 9 shows time domain filtering of multi-slice/multl-echo data using one-dimensional filter convolu- 
tions. 

Figures 1 0A-1 OD are views of exemplary reconstructions in the image space (Figure 1 0A) and time domain 
space (Figures 10B-10D), respectively, showing the effect of filtering in accordance with the method of the 
Invention. 

Throughout the drawings, like elements are referred to by like numerals. 



DETAILED DESCRIPTION OF THE INVENTION 



Figures 1 and 5 show an NMR phased array 10, such as described In Roemer et al. U.S. Patent No. 
20 4,825, 1 62, of a plurality of radio-frequency (RF) receiver coils 1 2 (coils 1 through Nc) defining an imaging volume 
for the NMR Imaging of a sample, such as for the NMR medical diagnostic imaging of a human spine. The sepa- 
rate surtece colls 1 2 are Identically configured and are arranged In doseiy-spaced relationship with overiapping 
fields-of-vlew (FOV), but with substantially no interaction between adjacent colls. The colls 12 are adapted as 
part of the NMR imaging process to simultaneously receive a different one of a plurality of NMR response slg- 
.■nalSAeachLexokedifnor ngggsg peiatedipjortioniofitheisampleAend 




35 Combination In the Image Domain 



As described in Roemer et al. U.S. Patent 4,871 ,969, tiij^giMmn^^ 
»theiindivJduaLcoiisjl^Nr-Jn^th.eir ^-r*— — ^» * — • * "^»^" 




45 



Assume that li(x,y) is the complex image obtained by reconstructing tn^aatanreceived from coil C|, and B| 
(x.y) is the RF magnetic field produced by coll Cj . The real part of B is the x component (in magnet coordinates 
as opposed to the screen coordinates of the image) of the transverse RF magnetic field and the imaginary part 
of B is the y component of the field. If noise correlations are ignored (which will have little effect on image quality) 
and all coils 12 have approximately the same noise, the combination of separate Images l| that optimizes the 
SNR in the composite image is given by 



50 



^ix,y) - ^ Iiix,y) Biix,y) 



(1) 



wher^B^^liiatbgggi ppj^^^ ^ .,.^^_-_^^.^^»^_p^_„__»„^^^ 

55 jTjhj^ comDlexiimageiisireallvithjeigi^^^ 

U x,v) = CS{x.ylBr(x,y) (2) 

Btie complex conjugate enters 
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equation (2) because increasing angles of th RF nnagnetic field are defined to be positive in the direction of 
rotation of the nuclei. Greater angles of the RF magnetic field con-espond to time delays (negative phase shifts) 
and thus the NMR signal Is proportiona l to the complex conj ugate of the RF magnetic field. 

ltlT§limaq (^g5a^^ 



lix.y) - X ^iix,y)Bi{x,y)\ - £ I Ji(x^y) II Btix,y) I 
li-i I i-i 



(3) 



Equation (3) gives a basicfonm usable in image reconstruction methods. Combining images using equation 
(3) is particularly convenient because the phase shifts of the individual receivers do not have to be known, and 
the image reconstruction programs do not have to carry the complex data. 

Figure 1 shows schematically the conventional process employed for combining the ^Q|^j"th® image 



domain. Thg|rjJMg|Si^^ 
ow^receiverMlS -"-'"^ ' 



30 



_ ^^^^ ^^^^^^^^ ^ SEJ^tillEroiinTiis^S 

f^||*5gg|{J|S}ggj^|^^ coil channel is separately subjected to 

transfonnation by processing means 19 and then combined point-by-point Into a single composite image at 
summation means 20 in accordance with equation (3). 

Combination In the Time Domain 

To derive the time domain filtering method of the present invention (Figure 5), it was recognized that the 
combined image obtained using the image domain method is simply the Fourier transform of the original tim 
dependent data. In accordance with equ ation (1), the optimal combination of images l (x,y) (i^^^that^^^^ 
SNR over the whole image) is obtained h.^gaffigigi^^ 
i^coilimaal^tiglfigliairoSg^^il^^ 

dSjTjgin^ for a single slic??™BR55iSBSffB^^ of the composite Image can 

S^lfSen by the two-dimensional convolution integral 
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A(trf to) « S( ( bL(tr.toWtr-Xr,t:o-to)clTodtr 



(4) 



i-1 — ^ 



whe relNgltslthMhirmbLeriofiCQiiS!vvi(t^ ^|il^^ 

■isitheiiSS^t^SiSiSj ^SSSr^^ is the inverse Fourier transform of the composite 

^SS^S^^S't^^^R^reaaourHni^ro^ac encode time t« . 

For a finite set of discrete samples, the inverse Fourier transform A(j,k) of the composite image is 
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A(j,K) = £ 1^ bi(l,m) vi(j-i,)c-m) 



(5) 



2 2 



where Vi(j,k) is a matrix of NMR voltages measured on coil i and b,0,k) is the discrete Fourier transfomi of the 
field map from coil i. The firs t and second argume nts are th e sample indices in the rea dout and ph ase encode 
directions, respectively, and ^^^^S^^^^^^^^S^^^^^^^S^^^^^^^^^^^^^ 



sgBectijSp im 

firstglance, equation (5) which combines the data in the time domain does not appear to offer any conri- 
putational advantages over equation (3) which combines the data in the image domain. According to equation 
(5), the total number of operations required to obtain a single k-space data point of the composite image is pro- 
portional to the number of pixels N in the image, where N=N<j> xNr. Thus, the number of operations required to 
construct the entire k-space representation of the composite image scales is a factor of N2 So, for large values 
of N, combining the data in im age space rather than in t ime space would appear to require far fewer compu- 
tations. This Is because tljggtaal5iTrQmieaGnklEjiilgh'aTin^^ etgaitoifSourieritpjisfo^ 
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<befd^conibining-^^ the number of operations for a Fast Fourier transform 

(FFT) seal s as a factor of N log(N) rather than N^. 

Th numb r of com^ations required for the convolution can be greatly reduc d, however, through the 
recognition that theRRfield map 21jshowTrinFi 
5 ancl^cafrtlius be suttablyTepresented in abbrevfiteia fonfn^ 

iabserved'that the inverse Fduri^^transfomi of tiie~fiell(i map is concentrated near the origin in t he time domai n 
(k-space) and thus the"l)li["1^) ten nsjn e guaBo^ be^tfuricatedto'a'kejt^ 

By way of exampleTF^ure 2 sRows ttie magnitude of a calculated sensitivity^profile of a typical s^gic^ 

Pcoil^The calculation is for a 40"arrFOy^itlTa~12l:if^^ coil located In a plane perpen-* 

10 dicular to the image. The main magnetic field is horizontal. The magnitude of its corresponding inverse Fourier 
transform is shown in Figure 3, which is a contour plot of the k-space representative of the field map of Figure 
2 . fiielronstant con tours are desjgnated-by-arbitrary-nu ^ 
maximum being affile prig lrTO nty the^enter 31x31~pixeIsof the magnitude are shown. Although thesensitivity\ 
^profile1(Figur e 2) qc<gjg^ 

15 wrflo mtud^ ^^ of j t^ filter fanction,profiig> 
[^iresponding to the-sensitlvlty^profile truncating the filter coefficient of Figure 3 by setting 

the magnitudes of the k-space representation^fTlgure 3 to zero outside the central 9x9 oixel matrix of points, 
placing i^Hammingwiricldw7(see, R.W. Hamming. Digital Filters prentice Hall] pp. 102-105) around the data 
to avoid ringing rand theh'Fourier t ransforming the res ult into image spacerA-visual compar te^^ 
20 r filter function profile of FigureTwi^^ difference, except^ 

^"hfe^rthe'colhwlres^th'fefn^^^ 
^tion—The error-near the wires (located at the intersection of the side lobes and the central region of senslth^ity) 
occurs t)ecause the RF magnetic field varies ra p|d_ly_there_a nd thus contains hJ gh^spatlaLfrequencies. A^jr^ 
^*=ofnttie coiirtRe"RF"fiera varies slowly and t he 9x9 filter kernel m atches the profjlewellr^ 



25 ^ — ^Using litjnrated'corivolutior^ to represent the field maps in k-space greatly reduces the number of 

operations necessary for effecting the conyoiutjon from one that s cales as a factor of to one that scales as 
a factor of NN|c, where Nk is the number of coefficients in the filter kerneL(e.g.. 81 for the 9x9 matrix above). 

Except for an arbitrary constsmtrthe'tiijhciated involution sca^ same as the Finr, 

Figure 5 shows a system and process in a ccordanc e with the invention fpfcombining the separate^itlda^ 

^ . _ _ — — _— — — - — ' r* ■ •..•■y,v ■---■v - . • -i^^i*i%l!<a^§^iPK^^§^3^JS,^fl!i^ 

30 \ from an NMR anray to~~dbtain a compositejimge with good SNR resQluUon in the time domainf using the above 




o 



„, , -, _ _„ ^ 

ro^MfijS^5Xt^ff5ht end of each coil channel 14' has its own receiver 15 and A/D circuitiy 16 for receiv^ 
ingand,digitlzingjh e^ separately received signaj. ^the same as for the corresponding channels 14 of the image 
domain processing set-up of Figure 1. But instead of storing the separate NMR images of each channel in a 
separate memory location 18 and Fourier transforming at 19 prior to summ ing , as done in the syg tem of Fi gur 
35 1 , the arrangementof FigureSrinaccordance witlvthe^invention. filtersthe_^data with a fieldjMp filter 23 as~ir7 
is-acquired,-and*tH^*sfims§*tli^^ 20' prior to stonng a pretransformation con>f 

binediimiai9'g"ir!>*IPW^gR^^ 

-means-25-to>S'ig||^^r^composite~image. Theliiters^S^rovide the weighting necessary for summrng th^ 
r separate contributions from the channiBlsjl41>& give^¥g6bli~SNRl^esolutiorr^^ image. The filters 23 

40 ^-"perfoffrTthe operations defirfed by equation (5). In contrast to the image domain data combination method 
employed by the system of Figure 1 , only o ne (or, possibry, two to obtain a uniform noise image) ff ourier trans— 
rfonmation-is-requirad arthe end'of the scanning op eration to produce thej sombjned.image. Thus the^^Sa blnainX 
^filtefin§t§ghlB?mfe^of FIgura the large time dela ys from end-of-scan to g g tjrnage^appearan cejnherent 

in the process-employed byThe syste m of Figure 1 . Moreover, th^ ^nsformation proce ss for the arrang^m ent^^}* 
45 shownirnggure Sjgndep^ number of ,c oils>utilized..^so, Irfcontraitto the Figure 1 image^omain 

approach, flie time domain method employ ed by the sys tem ofFj gure 5 of the invent ion_h as the additioria l 
advantage that the^dJata^i^^^^iiibli^ 
(^capaci tylieisessary for each coir dhahneir p> 

50 System Architectures 

Various system architectures are possible for implementation of the time domain image reconstruc tion 
method of the invention/lf the'k^space representaliorr(Figure 3)"of th~'fi~lcJ map for each coil 12 is reducecl to 
a kemet size of-N^x-N^in the readout and phase encode directions (s e Figure"4)rthe k-spa ge re presentatiop 
55 of the com posite imaoe-mav be Qiven^bv— — 
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A(j,}c) = ^ ]^ bi(l,m) Vi(j-l,k-m) 



(6) 
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According to equation (B^^acljudataiDon^^ 



ri3BSiteiktsuat>e|i^i ^ia uii , 

Figure 6 is a pictorial representation of truncated convolution for a two-dimensional single slice image. The 
digitized signal 28 from a single coil 1 2 has a number of data points 29, 30 of magnitude Vp , Vn+i . etc. , corre- 
sponding to samplings of the analog v oltage signal 28 taken at su ccessive time Inters/ats tn , t^+i, etc. jjjjgHatej 

iBS^sJfBRitributes to a rectangular subregion 32 shifted in the data memory 24 in the readout direction by 
one column. For each new phase encode step, the region moves down by one row. 

A straightforward hardware embodiment of filter 23 for implementation of equation (6) is depicted in Figure 
7. the arrangement uses discrete components and needs a modest amount of memory to temporarily save the 
most recent Nc* x data points from each slice, echo and coil. The multiplication element A performs the mul- 
tiplication in equation (6). The two Innenmost summations of equation (6) are done by summation element B 
and the te mporary register 36. The outemriost sum is perfonmed by summation elem ent 20jg|^ii^adij|^^^ 
inetfc1:4y ^^K§ mSSiaSs!Bi3i^^ 

WmmS r^lewSr!^ functions to save the data points 29, 30 for the successive time increments tn, tn+1. 
etc., until a sufficient amount of data is accumulated to complete one i,k point in ttie composite matrix. Assuming 
ttiat k-space is covered in a linear fashion, this occurs during tiie Nco 'th Phase encode step. As each subse- 
quent data point enters from the left (viz.. at the input end), the appropriate Ncr x Nc* data is exfracted from 
filter memory 34 and is multiplied by the corresponding filter coefficients stored in a convolution kernel memoiy 
35, summed by use of a temporary register 36, and then sent to an output port 37 as a weighted input to the 
composite summation means 20'. To avoid overflow of the input memory 34. a line of readout data must be 
completed and sent to the output port for each new line that entere. To accommodate multi-slice data, the mem- 
ory must save the most recent NcO x Nr data points from each slice and echo. 

One of the disadvantages of the architecture of Figure 7 is that each coil requires a significant amount of 
fast memory. The depicted memory 34 for each coil channel 14' has a size of H x Ncr x Ne x Ns x 4 bytes, wh re 
Nr is the number of readout points. N^ is the number of filter points in the phase encode direction, N© is the 
number of echoes, and Ng is the number of slices. The convolution kemel memory 35 has a memory size of 
Ncr X No* X Ns x 4 bytes, where No^, is the number of filter points in the phase encode direction. To keep up with 
the incoming data, the memory should operate at a rate that is N^r x Nc* faster than the entering data. Assuming 
a 100% duty cycle. 512 readout points, an 8 msec readout time and a 16x16 convolution kernel, the temporary 
storage memory data rates should be 16 MHz. A 32x32 point kernel would require a 64 MHz ""^JJig^BiBi^^p^ 
|g^[jgFging|the^^ 

EBg§gi^Oi^m^S rngure?siT^^ time domain data 

received fromeacif con 12 is operated on immediately as it entere from the left or input side of channel 14'. As 
each data point enters, it Is multiplied by a total of Ncr x Nc* differentfilter coefficients from tiie convolution kemel 
memory 35. The results of these multiplications are delivered on Inputs 37' to be summed at summation means 
20" over the coils 12 and then added to an output memory 34' inMreadHj^Q^ with the 

airangement shown in Figure 7, the data rates to and from tiie filter ^SlxS^^^SS^^lmmfmfnory 34' should 
be Ncr X Nc* times the incoming data rate. 

In the system of Figure 8, the output memory 34' is chosen to be a different element ttian that of the total 
composite memory 24 (Figures 5 and 6). In principle, ttiose memories could be constituted by the same physical 
memory, but splitting the memory into two sections may have cost advantages. A faster memory 34' at the out- 
put of summation means 20" is usually more expensive, but less of it is needed compared to tiie composite 
memory 24 (shown in Figure 5). 

The architectures shown in Figures 7 and 8 do not exploit any special purpose filt er chips t hat might be 

lakconyjolutionsweryi^ffiSientiSliaitd^ 
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available. A'nilimb'eV^fWbmtfantb'sWaK gTfchipsrthatieomPuteiOne^ 

hlteV^kigu^^ 

dalWaB^nBomSfflB^SfTBeTeSuce^ the phase 
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encode direction. For ach coil 12, a line of data in t he readout direc tion is passed through a one-dimensional 
convolution filter 38 a total of Nca» times5EactrpassIthrou^^ 

filter coe fftctents-obtained from th filter memory-and-corresponds to th -multiplication and^innem^^ 
maB^ofjS5!SS8f^^ th¥coils7ispToduced by summation 

5 means 20"'. Summation means 20'" also perfomns an addition for the read-mod ify-write operation which is 
equivalent to the second summation in equation (6). The signals at the outputs 37" of the filters are summed 
by summation means 20"'. and the results are added to the k-space storage memory 24 for the composite 
images. For each readout line that enters from the left or input side, Nc4> readout lines in the composite stage 
memor y 24 are nTOdifie d, Jhe same haitlwai^ca bychang- 

10 iiTigtheo^e^^ 
/ po^r^'memo^ 

One proposed digital signal processing chip employs 32 multipliers that operate in parallel and can con- 
volve a 16 bit-by-16 point complex kemel with a 16 bit data path at a rate of 800 ns per complex l/Q pair. For 
multi-slice data requiring a 16 point-by-16 point kemel and acquired at a 100% duty cyde, this corresponds to 

IS a maximum digitization rate of 12.8 ^sec. For conventional multi-slice imaging, 512 complex l/Q pairs can be 
typically acquired In 8 msec or 1 6 i^sec per point. Thus, it can be seen that a single filter chip may be sufficient 
for each coil channel 14', with some spare for overhead such as reloading the filter coefficients. 

For three-dimensional and high speed imaging, a number of chips would have to be placed in parallel to 
keep up with the data rates. An alternative approach for high speed imaging is to trade off the average rate 

20 with the instantaneous rate. Generally, high speed imaging uses small matrix sizes because of poor SNR. It 
may therefore be more cost effective to add a modest amount of memory at the input of the filters instead of 
coupling too many filter chips in parallel. The memory could then be filled at a high Instantaneous burst rate 
and emptied at a reduced average rate, giving the system time to filter the data. 

For three-dimensional imaging, the incoming data rates are the same as for multi-slice data, but another 

25 dimension of filtering is required. Since the data is acquired over many minutes with neariy 100% duty cycles, 
huge amounts of data are processed, it is therefore not practical to place temporary memory in front of each 
filter, and thus the data must be combined as fast as it enters. A three-dimensional image with 512 readout 
points and an 8 msec readout time would require 16 filter chips to keep up with the incoming data rate. A reduc- 
tion in the filter kernel from 16 to 8 in the two phase encoding directions changes the data rate through the filter 

30 by a factor of four and this only four chips would be required. This may cause some degradation in the SNR 
(initial indications are not much) but the SNR will still be better than one could obtain without using the phased 
array. 



Example 
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^ A four-cotlarrayv7as used to demonstrate the methods for sin gle slice sagflitai imaging of the human sptne.^ 
\The^arrayjvas.made-of~-1 2,cm jcoi|s_oyCTlapping in a row in a man ner simj iar^to that-shown-in-Figure-4 of-the'^ 
't62_patentJQiefour^ 12 were placed beneathlhe patiinfln a linear array running in the vertical direction. 
[^Each coil T27iadlts^vim^ image FOV was 40 cm with a composite matrix size of 

40 ^ 512x512 pixels. Afterthellafawas~sepa^ It was combined in the time domain with filter kemel 

sizes ranging from 1 x1 (a simple sum) to 9x9. For comparison, the data was also combined in the image domain. 

Following the general procedure described in the '162 patent, the data from each coil wa s conribined in the 
image domain_bv_separatei v Foll rl e^ from each coil 12, weighting each imagFby its.cor> 

^^sp onding field map im age,^nffByeffsum'm was then normalfeed into a uniform sct- 

45 sitivi ^ inrgge. Figure 10A show sgFresuitant 51 2x51 2 sagg ital image of the spine. 

[ To combine the data in thelirnellomairMJtilizing^ method'of the"present~invention,-the.filter-ooefficients 
were determined by caicufating tlie cbitiplex fi'^ia maip phase) for each coil over the^full 512x612/ 

image matrixTEx^i^ssiorTS^^ 

placedinrtlie centeriine of the coil. The complex field map for each coil was then inverse Fourier transformed 
50 into k-space and truncated to the desired size. To avoid ringing in the image, the resultant filter coefficients 
were windowed in two dimensions with a Hamming window of the type described in R.W. l-lammtng, . Digital 

Filters , supra. 

Tib avoid interference pattems m the imagerthe raw data wa^ GOTO different phases and^aihs 

of the NMR^^^ do thisv^i^tirah^f lotfl^oxiffi^ ^ced^bout..^r-^ 

55 .^xrrvabove ea^h suifI8§^lP^tfe~transmit loop wasdrivef^ the loc^roscillator used bylhe NMR receivers 
^^nd thiswas^hase^lockecJtb the receiver. The resultant signal amplitude and phase measured at the output 

of the receiv rs were used^to calibrate eac h channel. 

Using equation (6)fthefrun^tted1filter kemel yyas theacx)nvolved with the phas jDOfte^^ 
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r^resuite were then~summed bverthe coHs and^Fourier transformed.^^^tain a uniform.sensitivity imag^. the 
filter kernel was convolved with tfteoonjugate of the filter kern I and summed over the coils. This normalization 
map was also Fourier transfomned and divided pixei-by-pixel into the image. 

Figures 10B, 10C and 10D show the results of a simple sum, a 5x5 point and 9x9 point kernel. The simple 

5 sum Image (Figure 10A) corresponds to an image obtained using a single large coil. As expected, the simple 
sum image had poor SNR (approximately two times lower than the point-to-point weighted sum image of Figure 
10A) and poor suppression of motion and wraparound artifacts. The 5x5 kernel image showed significant 
Improvement The 5x5 Image (Figure 10C) was developed by passing the time dependent data from each coil 
12 through a two-dimensional filter with a 5x5 kernel before summing. A faint wraparound artifact is visible at 

10 the top. The 9x9 image (Figure 10D) was constructed similarly using a 9x9 kernel. The 9x9 result had almost 
the same quality as the composite image (Figure 10A) developed using the image domain techniques, except 
for minor differences near the edges of the image. 

The differences of the edges between the images obtained using time domain (Figures 10B and 10D) and 
Image domain (Figure 10A) methods are due to the wraparound of the filter. Ideally, the filter corresponding to 

IS the coil at the bottom of the image should have no significant contribution at the top of the image. However, a 
9x9 filter can be made to roll off In only about 1/9th of the Image FOV and. thus, some wraparound is unavoid- 
able. Combining the data in the image domain, however, allows one to filter to the nearest pixel or 1/512 of the 
image FOV. To obtain exactly the same result in the time domain would required a 512x512 filter kernel. 
In the above example, the filter coefficients were determined for the time domain filter hardware by cal- 

20 culating the RF magnetic field for each pixel in each slice for each coil. The results were then Fourier transfor- 
med, truncated, and then windowed. This method may not be fast enough jiowever, to be practi cal in the clinic a[ 
environment. ^Thefilter-coefficientsare^functlon^^ 
system-selectSfRi'SfSSSIS^^ 
C^Ttwo-after tfe Slectio^ be read yjgtake_data, y^g the method of the exampi , 

25 Vtodetermine the coefHcients for 30 si ices of 512x512 pixel images using thefour-coil anray. 120 two-dimensional 
512x512 complex inverse Fourier transforms would be required. If each transform takes 3 seconds on an anray 
processor, this part of the computation would take 6 minutes. The calculation of the magnetic fields would prob- 
ably add a few more minutes, thus creating a built-in delay of about 10 minutes, which might be unacceptable. 
A more rapid means of calculating the filter coefficients from the known positions or the RF receiving colls may 

30 be needed. 

There are a couple of approaches to speeding up the process of filter calculation. One approach is to pre- 
calculate and store filter coefficients for common slice locations, but this might be too restrictive. Anoth r 
approach involves a Fourier transform on a smaller grid. Since the resultant set of filter coefFicients in k-space 
will be truncated, it is not necessary to sample the RF magnetic field at each and every pixel in the image befor 
35 Inverse Fourier transformation. The matrix size for coefficient calculations can thus be reduced from, say, 512 
to perhaps 50 or fewer pixels. In the above 30 slice calculation, this would decrease the time by a factor of 1 00. 

r ^ J A furtberTJp)|:o3ch-in^ Sl^iP^^- 

I J/ (^<^ll',%JRFw!Xiagi2^^ position oMhe coil. t1|feTsl3o1^^^ 

I and Bz,_o f the ma qnetic_fieid.yEach component of the "RFlnagnetic~field"is"thenTnverse Fourier transform d, 

40 truncated, and saved on disk for later use. Since rotatio ns in real space are simple rotations in k- space and 
translations in real sp ace are phase j hift s in k-s pacefthe set of filter coefficients can be derived for any c<^l 
.location or orieritatiorTfrom this original stored set. Forjhree^ini^ 

'rotated and translated in k-space and therTwIrTd^owed to the desired size. For multi-slice data, the three-dimen- 
sional k-space data can be first rotated and then translated according to the slice position. TheUata can then 
45 ^bacollapsed into two dimfelTSiOns before windowing.-Such methods involv 
Cmatrice$.^Jhey-can-beaccomplished-quickly.. 



In yet another approach, the NMMR image data itself can be used as the basis for determination of the 
filter coefficientsJ]Tiisjsam and relies on the fact thatlhelmage 

f^fe^lf isTmeasureoTtl^^ 
50 '^"^^ oiii^liBif ^SSiS 

~There~has thus been described a method for combining the data from the separate coil channels of an NMR 
phased array in the time domain using filters to produce a composit Image having high SNR throughout the 
image. When compared with methods that combine the data in the image domain, substantial reductions in 
the reconstruction time and the amount of memory required in the NMR imaging system are realized. In this 
55 way, systems using many coils can be made more practical. 

Thos skilled in the art to which the invention relates will appreciate that other substitutions and modifi- 
cations can be made to the described embodiment without departing from the spirit and scope of the invention 
as described by the claims below. 
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Claims 

1 . A method for combining NMR response data of a samp! from a plurality of RF receiver coils of an NMR 
phased array in the time domain to form a composite NMR image, comprising the steps of: 

<a) receiving at each of the RF receiver coils a different one of a plurality of NMR response signals, 
each of the signals being evoked from a portion of the sample within a field of view of a respectiv one 
of the receiver coils; 

(b) conditioning each of the response signals to develop a plurality of data point signals corresponding 
to the magnitude of each of the respective response signals from each of the receiver coils at successiv 
time intervals; 

(c) convolving each of the data point signals by a time domain representation of a field map of the res- 
pective one of the receiver colls generating the corresponding one of the response signals; and 

(d) combining the signals obtained by the step of convolving on a time domain point-by-point basis to 
produce a time domain representation of the composite NMR image of the sample. 

2. The method of claim 1 wherein adjacent ones of the response coils have overlapping fields-of-view. 

3. The method of claim 2 wherein flie response coils have substantially zero interaction. 

4. The.method of claim 1 wherein the step of receiving includes the step of detecting each of the respons 
signals frem each respective one of the receiver coils substantially simultaneously. 

5. The method of claim 1 wherein the step of conditioning includes the step of digitizing each of the response 
signals to fonm a plurality of different series of discrete ones of the data point signals. 

6. The method of claim 5 wherein the step of convolving includes the step of reducing each time domain rep- 
resentation to a central kernel size. 

7. The method of claim 6 wherein the step of convolving is perfonmed in accordance with the relationship: 

8. The method of daim 1 wherein tiie step of convolving includes the steps of obtaining an initial NMR image 
representation of the sample from each of the RF receiver coils and substituting the initial NMR image 
representation for the time domain representation of each receiver coil. 

9. The method of claim 1 wherein the NMR composite image comprises a spectroscopic image of the sample. 

10. A method for combining NMR response data from receiver coils of an NMR phased array in the time domain 
using filtering to form a composite NMR image having good overall SNR resolution, comprising the steps 
of. 

(a) providing a plurality of dosely-spaced RF receiver coils, with adjacent coils having overlapping 
fields-of>view and substantially no interaction; 

(b) substantially simultaneously receiving at each one of the coils a different one of a plurality of NMR 
response signals, each evoked from a portion of the sample within the field-of-view of that coil 

(c) digitizing the response signal of each coil to provide a plurality of different series of discrete data 
point signals corresponding to the magnitudes of the response at successive time intervals for each 
coil; 

(d) temporarily storing an m x 1 number of said data point signals of each of the series of signals; 

(e) multiplying the temporarily stored signals by a con-esponding m x 1 numb r of signals corresponding 
to the digitized time domain representation of a field map of the coil associated with the series reduced 
to a cenbral m x 1 kernel size; and 

(f) combining the multiplied temporarily stored signals of each coO, on a time domain point-by-point 
basis, to produce a composite NMR image of the sample. 
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11. Apparatus for receiving and combining NMR respons data from a plurality of RF receiver coils of an NMR 
phased array to form a c mposite NMR image of a sample, comprising: 

a plurality of receiver circuits each connected to a respective one of the coils for rec iving a corre- 
sponding NMR response signal, each response signal being evoked from a portion of the sample within 
5 the field-of-view of the response one of the colls; 

an anaiog-to-digital converter connected to each receiver circuit for providing a series of discrete 
d ig Ital signals corresponding to a succession of time space data points representative of the response sig- 
nals from each one of the coils; 

means associated with each coil for providing an m x 1 number of digitized signals corresponding 
10 to the digitized time domain representation of a field map of the associated cofl reduced to a central m x 

1 kemel size; 

means connected to each converter and to said means associated with each coil for temporarily 
storing each last m x 1 number of said series of discrete digital signals; 

means for convolving each series of discrete digital signals with said digitized signals of said time 
15 domain representation of each field map to provide a convolved signal series; and 

means connected to said temporary storing means for combining said convolved series to provid 
a time domain representation of a composite NMR image of the sample. 
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@ Time domain filtering for NMR phased anray imaging. 

(g) A method and apparatus for combining NMR 
response data of a sample from a pturalify of 
dosely spaced RF receiver colls of an NMR 
phased array in the time domain to form a 
composite NMR image wherein each of the RF 
receiver coOs receives a different respective one 
of a plurality of NMR response signals, each of 
which is evolced from a portion of tiie sample 
within a field of view of a respective one of the 
receiver coils. The response signals are con- 
ditioned to develop a plurality of data point 
signsds corresponding to the magnitude of each 
of the respective response signals from each of 
the receiver coOs at suoces^ve time intervals. 
The data point signals are convoked by a time 
domain representation of a field map of the 
respective one of the receiver coils generating 
the corresponding one of the response signala 
The convoh^ed signals are combined on a time 
domain point-by-point t>asis to produce a time 
domain representation of the composite NMR 
image of the sample. 
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